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Notation 


ALPHA 

- 

see  a 

AXPS 

- 

intake  duct  static  pressures  (kPa) 

B 

- 

atmospheric  pressure  (mba) 

BETA 

- 

see  8 

DC-60 

- 

see  DCu, 

DOt,  ii 

- 

distortion  factor 

DP 

- 

see  AP 

EPS 

- 

see  e 

GDPTPO 

- 

array  containing  the  total  pressure  ratios  at  the  engine 
face  in  grid  format  as  shown  in  figure  11 

JPRESS(i) 

- 

array  containing  59  scanivalve  pressure  readings 

m 

- 

mass  flow  (lb/'sec) 

m 

- 

number  of  grid  points  in  radial  direction 

MDOT 

- 

see  th 

n 

- 

number  of  grid  points  in  circumferential  direction 

Pn 

- 

atmospheric  pressure  (kPa) 

P,„ 

- 

true  equivalent  pressure  from  scanivalve  (kPa) 

P. 

- 

static  pressure  (kPa) 

P, 

- 

engine  face  total  pressure  (kPa) 

P, 

- 

engine  face  dynamic  pressure  (kPa) 

P, 

- 

tunnel  freestream  total  pressure  (kPa)  (static  +  dynamic) 

Pi 

- 

pressure  before  venturi  throat  ( In.H<0 ) 

P  (’j) 

- 

PTM(t.j)  in  1-D  format 

PM  AX 

- 

maximum  engine  face  total  pressure  ratio 

PM  IN 

- 

minimum  engine  face  total  pressure  ratio 

PO 

- 

differential  tunnel  total  pressure  (kPa) 

PS 

- 

differential  engine  face  static  pressure  (kPa) 

PSA'' 

- 

average  of  six  static  pressure  ratios  around  engine  face 

PSPO 

- 

P./P„ 

PT 

- 

differential  engine  face  total  pressure  (kPa) 

PTAV 

- 

average  of  thirty  P,/P„  values  at  engine  face 

PTM(t,i) 

- 

interpolated  and  measured  engine  face  total  pressure  ratios 

PTMIN(U) 

- 

PTPO(i)  in  2-D  formal 

PTPO 

- 

see  PTAV 

PTPO(i) 

- 

thirty  engine  face  total  pressure  ratios 

Ptu 

- 

tunnel  freestream  total  pressure  (kPa)  (static  L  dynamic) 

PT60 

- 

minimum  of  (average  P,/  P,  over  60°  sectors  at  5°  intervals] 

PT60PS 

- 

sector  in  which  PT60  occurs 

PVAV 

- 

average  of  thirty  Px  / P,  values  at  engine  face 

PI 

- 

see  P\ 

R 

- 

current  value  of  annulus  number 

RO 

- 

SCG  px ,  || 

iii 


SECAVG 

T„ 

T, 

TA 

TH 

T1 

VEL 

A\, 

XMDOT 


sector  average  total  pressure  ratios 
atmospheric  temperature  (K) 
temperature  in  the  venturi  (°C) 
see  T„ 

current  value  of  circumferential  position 
see  T, 

tunnel  freestream  velocity  measured  bv  scanivalves  (m 's) 
scanivalve  pressure  reading 
mass  flow  parameter  m  \/T\,j  P, 


a 

& 

AP 

e 

A 

P,.„ 


pitch  angle  (deg)  -  nose  up  positive 
yaw  angle  (deg)  -  nose  to  port  positive 

pressure  difference  from  Pt  position  to  the  venturi  throat  ( In.ff.O ) 
expansibility  factor 
venturi  pressure  ratio 
venturi  air  density  (lb/ ft') 
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1  Introduction 


A  test  program  was  carried  out  in  the  ARL  low  speed  wind  tunnel  on  a  model  of  the 
Jindivik  (Figs  1.2).  to  evaluate  the  incorporation  of  an  auxiliary  intake  designed  to 
improve  performance  at  low  speeds.  A  FORTRAN  program  CALJIN  was  written 
on  the  DEC1  PDP  11/44  to  process  the  data  and  output  the  results. 

A  pressure  rake  containing  thirty  total  pressure  probes,  five  radial  positions  at  six 
circumferential  locations  (Fig.  3).  was  located  at  the  engine  face.  Pressure  mea- 
suiements  were  obtained  via  two  scanivalves  connected  to  the  pressure  rake. 

The  pressures  were  measured  using  the  wind  tunnel  scanivalve  system,  and  the 
output  data  was  acquired  using  the  standard  wind  tunnel  data  acquisition  program 
DATAIN  (Ref.  1).  Due  to  the  size  of  the  available  memory  it  was  not  possible  to 
install  CALJIN  as  a  subtask  of  DATAIN.  therefore  the  raw  data  were  output  to 
a  file  for  processing  by  CALJIN  independently.  Although  this  method  is  slightly 
cumbersome,  many  of  the  raw  data  files  needed  adjusting  because  of  faults  in  the 
automatic  data  acquisition  process.  These  faults  usually  involved  the  loss  of  one 
character  in  a  data  line  and  it  had  to  be  inserted  manually. 

The  results  were  expressed  as  pressure  ratios  and  they  were  output  in  graphical  and 
numerical  format.  The  test  results  were  presented  on  3  different  types  of  plots  as 
follows. 

1.  A  graph  of  the  ratio  of  intake  static  pressure  to  tunnel  freestream  total  pressure 
(static  4-  dynamic)  vs  axial  location  along  the  duct  away  from  (lie  engine  face 
towards  the  intake  lip.  (Fig  4) 

2.  A  graph  of  the  average  ratio  of  total  pressure  at  the  engine  face  to  tunnel 
freestream  total  pressure  ( Pi/P„ vs.  mass  flow  parameter  my/X,/ P„  was 
plotted  from  data  for  a  range  of  mass  flows.  (Fig.  5) 

3.  A  cross  sectional  total  pressure  ratio  contour  plot  at  the  engine  face,  as  shown 
in  figure  6. 

In  order  to  obtain  more  data  points,  (i.e.  more  than  the  30  readings  taken),  needed 
to  plot  contour  levels,  it  was  necessary  to  interpolate  values  between  the  measured 
pressure  points.  Initially,  linear  interpolation  methods  were  used,  but  the  contour 
plots  were  not  very  smooth  and  the  results  were  not  accurate  enough.  A  more  accu¬ 
rate  cubic  spline  interpolation  scheme  was  therefore  used  to  produce  the  interpolated 
pressures  needed  for  the  contour  plots. 

In  addition  the  numerical  data  were  presented  in  tabular  form  (Table  1).  The  plots 
were  output  on  the  IIP"’  7220T  plotter.  Although  the  plotter  was  slow  the  quality 
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and  clarity  of  the  plots  produced  was  good,  as  shown  in  the  examples  reproduced 
in  figures  4.  5,  and  6. 


2  Program  Structure 

The  program  structure  is  illustrated  in  figure  7  in  the  form  of  a  tree  diagram,  showing 
the  major  tasks  (programs)  and  subroutines.  The  source  code  for  these  files  can  be 
found  in  the  directory  DU0:f200.10j  of  the  PDP  11/44  (see  Appendix  A).  The  only 
exception  being  HPPLT.  which  represents  the  HP  plotting  library  routines,  which 
are  located  in  DU0:(1,1|HPPLTLIB. 

The  program  JINDIV,  which  is  essentially  a  menu-driver  to  invoke  each  separate 
task,  is  divided  into  four  main  sections: 

•  CALJIN  -  the  calculation  and  processing  of  the  measured  data 

•  PLTPS  -  the  static  pressure  graph. 

•  PLTPT  -  the  mass  flow  graph. 

•  PLTCNT  -  the  contour  pressure  plot. 

The  division  of  the  program  into  these  sections  allowed  each  task  to  be  developed 
and  debugged  individually  without  any  connection  to  the  other  tasks.  PLTPS  was 
incorporated  as  a  subtask  of  CALJIN  after  the  completion  of  the  programming, 
resulting  in  three  distinct  tasks,  CALJIN.  PLTCNT,  and  PLTPT. 

Initially,  it  was  assumed  that  the  entire  data  processing  and  output  routines  would  be 
invoked  as  a  subtask  of  the  standard  data  acquisition  program  DATAIN.  However, 
due  to  the  memory  limitations  of  the  PDP  11/44  this  was  not  possible,  and  it  was 
necessary  for  the  programs  to  be  written  as  separate  modules. 

The  existing  program  NORMALSPM  (a  subtask  running  under  DATAIN)  was  mod¬ 
ified  and  renamed  NORMALJIN  to  output  the  fixed  data  and  scanivalve  measure¬ 
ments  to  a  file  DL1:JINXXX.JCS \  which  could  be  read  by  CALJIN.  During  the 
tests,  faults  occurred  in  acquiring  a  preset  data  value  (Preset  6)  from  the  data 
serializer,  therefore  this  was  manually  input  bv  the  operator  at  runtime. 


2.1  Program  JINDIV 

The  program  JINDIV,  which  is  a  separate  major  task  on  the  PDP  11/44,  sets  up 
the  main  menu  for  the  processing  of  pressure  data,  and  output  of  graphs  and  plots. 

‘All  files  are  in  directory  (100,7), unless  otherwise  stated,  and  this  will  not  be  included  in  the 
file  description. 
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Three  separate  tasks  CALJIN.  PLTCNT.  PLTPT,  are  invoked  through  this  program, 
with  event  flags  52.  53  and  54  respectively.  At  the  completion  of  a  selected  task  the 
program  returns  to  the  main  menu. 


2.2  Task  CALJIN 

The  program  CALJIN  is  the  main  component  of  the  computer  software  written  for 
the  Jindivik  test  program.  All  of  the  data  processing  is  carried  out  in  CALJIN 
and  its  related  subroutines,  and  the  processed  data  is  output  to  files  for  printing  or 
further  analysis. 


2.2.1  Input 

Input  to  CALJIN  consists  of  two  items  of  information. 

•  First,  the  user  must  indicate  if  the  test  configuration  to  be  processed  includes 
one  of  the  range  of  auxiliary  air  intakes  tested.  If  an  auxiliary  intake  config¬ 
uration  is  present  static  pressure  hole  number  10  does  not  have  a  measured 
reading  (Fig.  8).  and  this  must  be  taken  into  account  to  avoid  it  being  plotted 
incorrectly  as  a  data  point.  A  set  value  of  P./P,,  =  0.55,  which  is  well  outside 
the  test  range  is  assigned  to  hole  number  10  when  an  auxiliary  air  intake  is 
included  in  the  tests. 

•  The  second  input  is  the  name  of  the  data  fiie  to  be  processed. 

The  standard  data  acquisition  program  DATAIN  asks  for  a  data  file  name. 
Fhe  user  selects  a  file  name  in  the  form  JINXXX  where  XXX  is  the  sequential 
test  number.  The  subroutine  NORMALJIN  (modified  from  NORMALSPM) 
outputs  the  pressures  and  other  data  to  a  file  DLkJINXXX.JCS  which  is  then 
input  to  CALJIN.  The  user  only  needs  to  type  JINXXX  without  the  necessary 
extensions. 

The  format  of  the  file  DLkJINXXX.JCS  (Table  2)  is  as  follows: 

1.  The  first  6  lines  contain  the  59  seanivalve  readings  (two  18  port  scamvahes 
were  used)  in  integer  free  format,  followed  by  a  blank  line. 

2.  The  next  2  lines  of  data  are  the  7  Preset  values  read  in  from  the  data  serializer 
in  real  free  format. 

3.  The  next  5  lines  are  test  description  lines  in  character  format. 
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4.  The  following  7  lines  contain  the  title,  customer  name,  test  point  number,  data 
file  name  (DLlrJINXXX.DAT1).  date  of  file  creation,  configuration  file  name 
and  configuration  value'. 


2.2.2  Scanivalve  Pressure  Calibrations  and  Pressure  Readings 

Pressure  tubes  from  the  Jindivik  model  were  connected  to  two  SCANCO  48D  type 
scanivalves  located  under  the  working  section  of  the  wind  tunnel.  STATU  AM  ±5 
psi  differential  pressure  transducers  were  used  in  the  scanivalves.  The  transducers 
were  not  linear  over  the  entire  pressure  range. 

Calibration  of  the  transducers  revealed  a  linear  relationship  over  the  range 
-24 kPa  <  P,„  <  -5-24 kPa.  however  a  definite  increase  in  slope  is  evident  out¬ 
side  this  range  (Figs  9.10).  Linear  equations  for  each  of  these  regions  relating  the 
scanivalve  pressure  reading  A\, .  and  the  true  equivalent  pressure  Pm.  in  kiloPascals. 
were  derived  leading  to  four  equations,  two  for  each  of  the  transducers. 

The  calibration  equations  for  Scanivalve  D  with  Transducer  54820  are: 


P,„  =  0.0451 13A\,  -i-  10.692  for  .V.,  £  -769  (1) 

P,„  -  0.0451 13A\,  -  10.692  for  A'.,  >  769  (2) 

The  calibration  equations  for  Scanivalve  E  with  Transducer  57045  are: 

P,„  -=  0.040956.W,  4-  5.835  for  A'.,  £  -730  (3) 

P,„  -  0  040956 A'.,  -  5.835  for  A\,  >  730  (  I) 


A  known  pressure  of  21.0 kPa  was  applied  to  a  port  on  each  of  the  scanivalves  during 
the  tests.  A  calibration  factor  for  each  of  the  scanivalves  can  then  be  calculated  for 
the  range  of  values  -769  <  A-.,  <  4-769  (scanivalve  D),  and  -730  <  A\,  £  ‘-730 
(scanivalve  E),  and  used  to  determine  pressures  within  these  ranges. 

Scanivalve  pressure  readings  are  read  into  variable  JPRESS.  an  integer  array  with 
59  values.  Pressures  are  grouped  in  JPRESS(i)  as  shown  in  Table  3. 


Table  3:  Variable  JPRESS(i) 


i 

Description 

1.2 

3.4 

5,6 

7,8 

9. . ...38 

39.. ... 44 

45..  . .  ,59 

Atmospheric  Pressure 

Calibration  Pressure  of  21.0 kPa 
Tunnel  freestream  Static  Pressure 
Tunnel  freestream  Total  Pressure 
Engine  face  Total  Pressure 

Engine  face  Static  Pressure 

Intake  duct  Static  Pressure 

'This  is  the  raw  data  file  created  by  DATAIN 

’These  values  are  not  all  required  by  CALJ1N  even  though  they  are  written  to  the  data  file 
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The  pressures  are  output  as  ratios  to  the  tunnel  freestream  total  pressure,  i.e.  P,/ P„ 
and  PJP„.  The  differential  pressure  transducers  measured  the  difference  between 
the  pressure  on  a  given  model  pressure  tapping  and  atmospheric  pressure  P„  which 
are: 

PT  -  P,  -  P„ 

PS  =  P-  P„ 

PO  =  P„  -  P„ 

To  calculate  the  correct  pressure  ratios,  P„  must  be  added 
and  the  denominator,  so  that  the  total  pressure  ratio  is: 

P  =  PT  +  P„ 

P,  PO  +  P„ 

and  the  static  pressure  ratio  is:*' 

P  PS  +  P„ 

P„  “  PO  +  P„ 


to  both  the  numerator 

(5) 

(6) 


The  six  measured  static  pressures  around  the  engine  face  are  averaged  to  give 
(P,/P,, The  dynamic  pressure  ratios  can  then  be  calculated  using  equation  7. 


(7) 


2.2.3  Cubic  Spline  Interpolation  for  P,/Ptl 

Total  pressure  measurements  were  recorded  at  thirty  positions  at  the  engine  face, 
over  five  annuli,  each  annulus  containing  six  pressure  probes  (Fig.  3).  To  carry 
out  the  necessary  calculations,  and  to  produce  a  reasonable  pressure  contour  plot, 
intermediate  interpolated  pressure  values  at  5°  intervals  around  each  annulus  are 
required.  A  linear  interpolation  method  was  originally  implemented,  however  the 
resulting  contour  plots  were  not  very  smooth  and  the  correlation  with  sample  data 
from  the  static  tests  was  not  satisfactory.  A  cubic  spline  interpolation  was  used 
which  produced  more  accurate  results. 

The  cubic  spline  routines  are  found  in  DU0:[200,10|CUBSPL.FTN,  and  are  called  In 
both  CALJIN  and  CNTJIN.  The  two  subroutines  SPLINE  and  SPLINT  contained 
in  CUBSPL  appear  in  their  original  C  format  given  in  Reference  2.  along  with  a 
detailed  explanation  of  the  algorithms. 

Interpolated  and  measured  pressure  ratios  at  the  engine  face  are  stored  in  a  two 

dimensional  array  GDPTPO(R,TH)  where  R  -  1, _ 5,  TH  =  1,. ..,72  (Fig.  11).  A 

polar  grid  was  chosen  because  it  allowed  the  results  for  each  annulus,  or  for  a  given 
sector,  to  be  manipulated  easily,  as  described  in  the  following  section. 

''In  CALJIN  the  static  pressure  variables  arc  PS  and  AXPS. 
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2.2.4  Engine  Face  and  Duct  Calculations 


Processing  of  the  measured  data  and  the  calculation  of  parameters  such  as  the  dis¬ 
tortion  factor  £)(?(,<>,  mass  flow  rate  m  are  found  in  the  following  section  of  CALJIN. 

Distortion  Factor 


The  distortion  factor  is  a  measurement  of  the  variation  of  total  pressure  across  the 
engine  face,  and  is  defined  by  equation  8. 


DC,,, 


(8) 


where: 


(£)„„is  the 
(&)«,  is  the 


average  of  the  thirty  total  pressure  values  at  the  engine  face: 
minimum  of  (average  (P,/ P„)  over  60°  sectors  at  5°  intervals); 
average  of  the  thirty  dynamic  pressure  values  .u  the  engine  face. 


In  order  to  calculate  (P,/ P„)w„  60°  sectors  are  chosen  by  taking  the  GDPTPO(R.TH) 
values  for  R=l,.  ..,5,  TH  =  1  (Lower  limit).. .  .,TH=13  (Upper  limit)  (Fig.  11),  cal¬ 
culating  the  average  of  these  total  pressure  ratios  and  storing  the  solution  in  the 
variable  SECAVG(l).  The  lower  limit  and  upper  limit  of  TH  are  incremented  by  1 
and  the  average  pressure  ratio  found  again  (SECAVG(2)).  This  process  is  repeated 
until  all  72  values  of  the  sector  averages  are  calculated.  The  only  tricky  section  of 
this  algorithm  is  to  remember  when  wrap  around  of  the  sectors  occur,  i.e.  when  the 
upper  limit  of  TH  equals  73.  values  of  TH  equal  to  1  must  be  used,  and  so  on. 

The  minimum  of  the  72  sector  average  values  is  chosen  and  this  is  assigned  to  the 
variable  PT60  (  (Pf/Pi)(,n),  allowing  the  distortion  factor  DCc,  to  be  calculated 
according  to  equation  8. 

Mass  Flow  Rate 

Air  mass  flow  through  the  engine  is  calculated  from  the  venturi  readings  of: 

P i  -  pressure  before  the  throat  (ln.H,0)\ 

A P  -  pressure  difference  from  P,  position  to  the  throat  (In. ll>0): 

7’|  -  temperature  in  the  venturi  (°C). 


The  following  equations,  derived  from  equations  found  in  Reference  3,  are  used  to 
calculate  the  mass  flow  m  and  the  mass  flow  parameter  ;hv/7'„/P„. 


Venturi  air  density  (lb/ ft'): 

p,.„  =  0.73291 


0.0296  D  -  0.0737  P,' 
273.15  +  7) 


(9) 
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Expansibility  factor: 


c 


0.9517 


1  -  0.0483  A1  '-*• 


(10) 


where  the  venturi  pressure  ratio  is  defined  as: 


A  = 


1  - 


AP 

P,  -  P, 


Mass  flow  rate  ( Ibjsec ): 

m  —  3.1735  e  \jp, , „  AP 
Therefore  the  mass  flow  parameter'. 


P„ 


\/273.15  +  T, 

g  I  I  (>y)(> 


00 


(12) 


(13) 


The  values  of  (P,/P,), ,,,,  and  th\/T„/ P„  for  each  test  (i.e.  one  mass  flow  setting)  are 
appended  to  a  file  DU0:XMDYPT.DAT.  This  file  is  accessed  by  the  task  PLTPT 
when  a  series  of  tests  have  been  completed,  and  a  summary  of  the  results  is  plotted. 


2.2.5  Subtask  PLTPS 

PLTPS  is  invoked  from  within  CALJ1N  and  produces  the  Static  Pressure  ratio 
(PJ Pi)  vs.  Axial  Location  (from  engine  face)  graph.  The  static  pressures  are  passed 
from  CAIJIN  via  the  file  DU0:PSTAT.DAT  containing  the  test  number,  the  sixteen 
values  of  P«/P„  from  the  engine  face  to  the  intake  lip,  and  the  test  description. 

Standard  HP  plotting  library  routines  DU0:(1,1]HPPLTLIB  to  control  the  1IP  7220T 
plotter  are  used  to  produce  the  plot  file  DU0:PSGXXX.PLT,  which  is  immediately 
queued  to  the  plotter. 


2.2.6  Total  Pressure  Ratio  Contour  Plotting 

The  grid  numbering  system  for  the  engine  face,  used  for  t he  manipulation  of  tin- 
pressure  ratios  for  contouring  is  different  to  the  numbering  system  previously  ex¬ 
plained  in  section  2.2.4.  A  diagram  illustrating  the  variation  in  the  grid  orientation 
is  shown  in  figure  11.  The  pressure  ratios  are  converted  from  their  original  one 
dimensional  array  variable  PTPO(i)  i  -  1 . 30,  to  the  corresponding  two  dimen¬ 
sional  array  format  PTM1N(1,J)  1  =  1, _ 7,  J  =  1,...,7.  Pressure  data  is  then  in 

the  correct  format  for  input  to  the  subroutine  DU0:(200,10]CNTJIN. 

'The  units  of  P„  in  this  formula  arc  psi  not  kPa 
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In  the  radial  direction  five  pressure  measurements  are  taken.  The  contour  program 
requires  values  at  both  boundaries  (i.e.  the  inner  and  outer  circles),  and  these  are 
calculated  using  linear  extrapolation. 

Two  subroutines  are  contained  within  CNT.JIN.  one  called  CNTJIN  and  the  other 
INTERP.  INTERP  performs  the  cubic  spline  interpolation  on  the  pressure  data  by 
calling  CUBSPL,  and  the  results  are  stored  in  the  format  necessary  for  input  to  the 
contour  routines. 

After  the  interpolation  is  carried  out  on  the  measured  total  pressure  ratios,  the  com¬ 
plete  total  pressure  data  is  then  transformed  from  the  two  dimensional  coordinate 
format  PTM(i,j)  i  =  1,...,7,  j  =  1, . . . ,  73,  (Fig.  11)  to  a  one  dimensional  array 
P(tj)  according  to  equation  14. 


=  i  +  t  =  1,  —  ,7  j  =  l,..., 73  (14) 

where  m  =  number  of  grid  points  in  the  i  direction. 

In  order  to  plot  contours  it  is  necessary  to  choose  a  range  of  contour  values  for  which 
the  pressure  ratios  will  be  plotted.  As  the  total  pressure  ratios  for  all  tests  always  lie 
within  the  range  0.84  <  P,/P, ,  <  0.99,  the  method  chosen  was  to  find  the  minimum 
and  maximum  pressure  ratios  (PMIN  and  PMAX),  and  to  divide  this  range  into 
increments  of  0.01. 

A  file,  DU0:PINXXX.DAT  (Table  4),  is  created  to  store  the  interpolated  pressure 
data  and  the  required  constants.  The  format  of  this  file  is  as  follows: 

1.  First  line  contains  a,  /3,  tunnel  velocity  (m/s),  m\fT\,/P, ,,  (P,/Pn ),„,,. 
( Pi / P(i)(,n,  and  DC(, n. 

2.  Second  line,  m  -  number  of  radial  grid  points,  n  -  number  of  circumferential 
grid  points,  and  the  number  of  contour  levels. 

3.  Third  line,  the  actual  contour  level  values. 

4.  Remainder  of  the  file  contains  a  table  of  pressure  ratios,  one  line  for  each 
circumferential  position. 

2.2.7  Output  of  Results 

Numerical  results  are  output  to  both  the  screen  and  a  file  DUOrOUTXXX.DAT  in 
the  tabulated  format  shown  in  Table  1.  Two  copies  of  this  file  are  placed  in  the  line 
printer  queue.  The  only  difference  between  the  screen  and  file  output,  is  that  the 
file  output  contains  a  title  page,  detailing  the  test  description,  name  of  customer, 
configuration  and  data  file  names,  and  the  creation  date  of  the  data  file. 
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0 .  °1 Co£69 
0 .  =: : 03: 1 
C.  =112063 
C .  =  *.  1 25  £2 
C .  =1 122  =  4 
0.  =11 1644 
C. =111072 
0 . =1 1 1 022 
0. =11 1=33 
0.=114244 
O.=1l£2oo 
0.012*4637 
0.0122614 
0.0144836 
0. =157236 
0.=:->3221 
0. =183445 
0.0195905 
0. =207034 
0 .  =21 ©156 
0. =222594 
0 .  c22£©7 - 
0.=224719 
0 .  =21 °053 
0. =208212 
0.  =1=2584 
0.  =172764 
0 .  =1 -°3S3 
0.  =  122  =  47 
0 .  =094146 
0.00625-7 
0.0221749 
0.5=09348 
0 . 8=6 t=46 

0  .  £c2=.l  38 
0. 6904E22 
0.6675610 
0.6646622 
0 . 8622649 
0. £80051 4 
0.6780221 
0. £762018 
0.8746290 
0.8722* ©3 
0.8722752 
0.871=171 
0.8710529 
0.8708969 


TABLE  4  INPUT  FILE  FOR  PLTCNT 
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2.3 


Task  PLTCNT 


This  task,  which  produces  the  total  pressure  contour  plot,  is  invoked  by  choosing 
the  second  option  in  the  main  menu  of  JINDIV.  In  order  for  a  contour  plot  (Fig. 
6)  to  be  created,  the  raw  data  for  that  test  must  be  processed  first,  i.e.  the  first 
option  (task  CALJIN)  of  JINDIV  must  be  carried  out  before  the  second  option  (task 
PLTCNT). 

2.3.1  Input 

PLTCNT  requests  the  test  number  and  a  one  line  test  description  of  the  test  for 
which  the  user  wants  to  produce  a  contour  plot.  The  file  DUO:PINXXX.DAT  con¬ 
taining  the  pressure  ratio  data  and  constants  is  read  into  the  appropriate  variables. 
The  format  of  the  file  PINXXX.DAT  is  described  in  section  2.2.6. 

2.3.2  Contour  Program 

A  standard  contour  plotting  package  DU0:(200,10]CONTOUR.FTN  is  used  to  create 
the  pressure  contours  using  the  interpolated  pressure  data  from  PINXXX.DAT. 
Unfortunately  this  contour  package  does  not  label  the  contour  levels  and  they  must 
be  labelled  manually  as  shown  in  figure  6.  To  make  this  process  easier,  the  program 
CONTOUR  was  modified  to  alternate  the  line  type  of  the  contour  when  a  change 
in  contour  level  occurs. 

The  plot  commands  are  written  to  a  file  DUO:CNTXXX.PLT  as  they  are  calculated. 
After  the  program  CONTOUR  has  been  completed,  the  title  and  test  information 
is  appended  to  the  plot  file,  followed  by  the  plot  commands  to  draw  the  inner  and 
outer  circles.  In  addition  to  the  constant  pressure  contours,  dashed  lines  are  placed 
at  60°  intervals  to  show  the  positions  of  the  probe  rakes  around  the  engine  face. 

2.3.3  Contour  Plot  Output 

The  file  containing  the  IIP  plot  commands,  CNTXXX.PLT.  is  then  placed  in  the 
queue  for  the  HP  7220T  plotter.  Contour  plotting  is  the  most  time  consuming  part 
of  the  processing  procedure,  as  some  of  the  plots  may  include  up  to  fourteen  ronioui 
levels.  The  speed  of  the  plotting  could  be  increased  by  using  a  faster  plottei  or  a 
laser  printer/plotter. 


2.4  Task  PLTPT 

This  task  is  invoked  by  choosing  the  third  option  in  the  main  menu  of  JINDIV. 
PLTPT  produces  a  plot  of  (P,/P„)(MI>  vs.  m^%/ Pn  for  a  range  of  mass  (lows  (Fig. 
5).  After  a  series  of  tests  are  completed  PLTPT  is  run  '  summarise  the  output 
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of  the  tests.  Any  discontinuity  in  this  graph  most  likely  indicates  the  presence 
of  erroneous  data  points,  and  it  was  found  that  this  was  the  easiest  method  to 
determine  if  incorrect  data  was  input. 


2.4.1  Input 

The  test  number  and  test  description  are  required  as  input  to  PLTPT.  The  test 
number  chosen  is  the  corresponding  test  number  of  the  last  test  in  that  series  of 
tests.  The  values  of  (P, / P, and  m\/T\,/ P„  from  file  DUO:XMDYPT.DAT  (which 
is  updated  through  CALJIN),  are  read  into  two  arrays  for  plotting. 


2.4.2  Plot  Description 

In  addition  to  the  graph  of  (Pi/Pi\),„,,  vs.  m\/%,jP„  being  plotted  from  the  two 
arrays,  constant  RPM  lines  representing  the  12,000  and  13,800  RPM  engine  settings 
are  plotted. 

The  equations  of  these  two  lines  are  given  by  equations  15  and  16  respectively. 


m\fP, 

P, 


lily/T, 


(15) 

(16) 


2.4.3  Plot  Output 

The  HP  plot  commands  are  written  to  a  file  DU0:PTGXXX.PLT,  and  this  is  placed 
in  the  queue  for  the  HP  7220T  plotter.  A  copy  of  the  file  DU0:XMDYPT.DAT  is 
made,  and  the  new  file  name  is  DU0:PTGXXX.DAT. 

It  is  very  important  to  note  that  after  a  series  of  tests  are  completed  and  the  plot 
file  PTGXXX.PLT  is  output  correctly,  the  file  contents  of  XMDYPT.DAT  must 
be  cleared  manually,  prior  to  processing  additional  tests.  'Phis  can  be  achieved  bv 
edit  ting  the  file,  deleting  the  two  columns  of  data,  and  saving  a  new  version  of 
XMDYPT.DAT.  If  this  is  not  performed,  data  points  will  be  appended  to  the  end 
of  the  file,  and  future  plots  will  include  more  than  one  series  of  tests. 


2.5  Subroutines 

Two  additional  subroutines  are  not  covered  by  the  previous  sections  of  this  docu¬ 
ment,  these  are  CATJIN  and  INFILE. 
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CATJIN  is  called  by  each  of  the  tasks  in  order  to  set  up  the  name  of  the  data  file  and 
plot  file  based  on  the  test  number.  Depending  on  which  task  calls  this  subroutine  the 
file  name  returned  is  a  concatenation  of  the  required  characters  and  test  number. 
If  CATJIN  is  called  by  PLTPS  with  parameters  flag  =  3  and  test  number  —  23. 
CATJIN  will  return  the  file  name  PSG023.PLT. 

INFILE  (a  subroutine  of  DATAIN),  called  by  CALJIN,  reads  in  the  file  name 
(JINXXX)  and  returns  the  full  file  name  description  (DL1:JINXXX.JCS). 


3  Software  Installation 

The  task  DU0:(200.l0]COMJIN.TSK  must  be  installed  in  order  for  DATAIN  to 
recognise  the  task  name  JIN.  COMJIN,  a  subtask  of  DATAIN.  performs  the  acqui¬ 
sition  of  the  scanivalvc  pressure  measurements  and  outputs  the  values  in  the  format 
described  in  section  2.2.1.  The  only  differences  between  COMJIN  and  C'OMSPM 
(the  original  subtask  for  scanivalvc  pressure  measurement)  are  the  modifications 
made  to  the  subroutine  NORMALJIN. 

The  programs  JINDIV.  CALJIN.  PLTPS,  PLTCNT  and  PLTPT  must  be  installed 
from  directory  DU0:[200,10).  A  command  file  DU0:|200, 10jJININS.CMD,  which 
includes  the  commands  to  compile,  link  and  install  the  tasks,  can  be  run  to  carry 
out  the  correct  installation  of  these  programs  and  subroutines. 


3.1  Procedure  to  Install  and  Run  Software 

The  following  DCL*  commands  for  the  DEC  PDP  11/44  will  install  and  run  the 
program.  Bold  faced  text  must  be  input  via  the  keyboard. 

>SET  DEF  DU0:[200,10]  set  directory  as  DU0:[200,10] 

>@COMJIN  if  changes  have  been  made  to  NORMALJIN 

>INS  COMJIN  install  task  COMJIN 

>@JININS  install  tasks  JIND1Y.CALJIN, PLTPS, PLTCNT. PL!  P  I 

>SET  DEF  DU0:(100,7]  or  login  to  directory  DU0:[100,7| 

"♦♦Ready  to  run  DATAIN M " 

>R  DATAIN 

l:  Job  name:  JIN 


Digital  Command  Language 


x  -  configuration  no. 
Preset(6)  not  functional 

****Exit  DATAIN**** 

Set  up  the  HP  7220T  plotter 

****  Ready  to  run  JINDIV"** 


2:  Data  file  name:  JINXXX 
3:  Configuration  file  name:  CONFJINx 
4:  Acquire  pressure  measurements  if  online 
5:  Enter  Patm:  P„ 

6:  Output  printed:  N 


JINDIV 

1:  Option:  1 
2:  Hit  <return> 
3:  Option:  2 
4:  Repeat  1  — *  3 
5:  Option:  3 
6:  Option:  Q 


process  raw  data 

until  main  menu  is  displayed 

output  contour  plot 
until  a  series  of  tests  are  processed 

output  mass  flow  graph 
Quit  task 


After  quitting  from  the  program  JINDIV  and  ascertaining  that  the  data  is  correct, 
clear  the  contents  of  the  file  DUO:XMDYPT.DAT  before  running  JINDIV  again. 


4  Conclusion 

Computer  software  developed  for  the  low  speed  wind  tunnel  PDP  11/44  computer 
and  described  in  this  document  made  the  post  processing  of  data  from  the  Jindivik 
tests  both  quick  and  easy  to  perform.  A  range  of  mass  flows  were  set  and  the  data 
acquired,  and  within  twenty  minutes  (depending  on  the  data  being  ’error-free’)  the 
complete  results  including  pressure  contour  plots,  graphs  and  tabulated  data  for  a 
series  of  tests  were  ready  to  be  analysed. 

The  scanivalves  are  now  being  phased  out  and  a  new  PSI  pressure  scanning  system  is 
being  installed  which  will  enable  the  pressures  to  be  acquired  at  a  much  faster  rate, 
and  ’error-free’  on  the  first  acquisition.  This  factor  and  the  possible  replacement  of 
the  HP  7220T  plotter  for  a  faster  plotter,  will  lead  to  a  considerable  reduction  in 
the  processing  time  of  wind  tunnel  test  results. 
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Appendix  —  Location  of  Files 

Source  files  in  directory  DU0:[200,10j: 


JINDIV.FTN 

CALJIN.FTN 

PLTPS.FTN 

PLTCNT.FTN 

PLTPT.FTN 

CNTJ1N.FTN 

CUBSPL.FTN 

CATJIN.FTN 

CONTOUR.  FTN 

NORMALJIN.FTN 


Data  files  and  plot  files  in  directory  DU0:(  100,7]  during  testing: 

OUTXXX.DAT 

PINXXX.DAT 

XMDYPT.DAT 

PTGXXX.DAT 

PSGXXX.PLT 

CNTXXX.PLT 

PTGXXX.PLT 


Data  files  in  directory  DLl:(100,7]: 

JINXXX.DAT  -  raw  data  file  from  DATAIN 
JINXXX.JCS  -  output  file  from  NORMALJIN 


Copy  of  data  files  and  plot  files  in  directory  DL0:[  1 00.7]: 

OUTXXX.DAT 

PTGXXX.DAT 

PSGXXX.PLT 

CNTXXX.PLT 

PTGXXX.PLT 
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FIG.  2  EXPERIMENTAL  SETUP  OF  JINDIVIIC  MODEL 
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FIG.  4  GRAPH  OF  Ps/P0  vs.  AXIAL  LOCATION  (mm) 


TEST  NO.  :  293. 

VEL  =32.9  XMDOT  =  3.122 

ALPHA  =0.0  3ETA  =  0.0 

PTAV  =  0.954  PT60  =  0.909  DC60  =  0.3 

PROFILED  LIP  -  MASS  FLOW  5-30  M/S 

FIG.  6  P,/P0  PRESSURE  CONTOUR  PLOT  AT  ENGINE  FACE 
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Scanivalve  Reading  X  sv 


PIG.  9  GRAPH  OF  SCANIVALVE  PRESSURE  Pm  vs. 
SCANIVALVE  READING  FOR 
SCANIVALVE  D  WITH  TRANSDUCER  54820 
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FIG.  10  GRAPH  OF  SCANIVALVE  PRESSURE  Pm  vs. 
SCANIVALVE  READING  X„  FOR 
SCANIVALVE  E  WITH  TRANSDUCER  57045 
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